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ABSTRACT 

The first stars in the universe, forming at redshifts z> 15 in minihalos with masses of order 10^~^Mq, may 
leave behind black holes as their remnants. These objects could conceivably serve as "seeds" for much larger 
black holes observed at redshifts z < 7. We study the growth of the remnant black holes through accretion in- 
cluding for the first time the emitted accretion radiation with adaptive mesh refinement cosmological radiation- 
hydrodynamical simulations. The effects of photo-ionization and heating dramatically affect the accretion flow 
from large scales, resulting in negligible mass growth of the black hole. We compare cases with the accretion 
luminosity included and neglected to show that the accretion radiation drastically changes the environment 
within 100 pc of the black hole, where gas temperatures are increased by an order of magnitude. The gas den- 
sities are reduced and further star formation in the same mini-halo prevented for the two hundred million years 
of evolution we followed. These calculations show that even without the radiative feedback included most 
seed black holes do not gain mass as efficiently as has been hoped for in previous theories, implying that black 
hole remnants of Pop III stars that formed in minihalos are not likely to be the origin of miniquasars. Most 
importantly, however, these calculations demonstrate that if early stellar mass black holes are indeed accreting 
close to the Bondi-Hoyle rate with ten percent efficiency they have a dramatic local effect in regulating star 
formation in the first galaxies. 

Subject headings: cosmology: theory — galaxies: formation — black holes: formation 



1. INTRODUCTION 

The first stars in the universe may have left a population 
of black holes as their remnants. This possibility demands 
we understand how radiative feedback would have affected 
their environment and regulated accretion. Luckily, there is a 
relative amount of theoretical consensus regarding the nature 
of first star formation, offering a solid foundation for further 
investigation into their black hole remnants. 

The first stars probably formed at redshifts z > 15 in mini- 
halos with mass lO^' ^Mf;^, and were hkely massive, with 
30 < M,/Mq < 300 (lAbel et al.112001 iBromm et al.ll2002h . 
They influenced subsequent structure formation through 
their UV radiatio n heating the surro u nding gas and initiat- 
mg reionization (IWhalen et a l.l l2004t iK itavama et al.l 120041: 
Alvarez et al. 2006; Johnson & Bromm 2007; Abel et al. 
2007i;,Yoshida et al.,.2007MWise & Abelli2007M200 8). If these 



population III (Pop III) stars were less massive than 140 M© 
or more massive than 260 Mq, they likely wo uld collapse 
to for m a black hole at the end of their lifetime ( iHeger et al.l 
l2003h . 

The presence of these black holes could have had dramatic 
consequences. If a lOO-M© black hole is radiating at the Ed- 
dington limit with ten per cent radiative efficiency, beginning 
at z > 20, it will attain a final mass at z = 6.4 of greater than 
lO^M©, implying that accretion onto these black holes is a 
viable explanation for the z^6 quasars (e.g.. lHaiman & Loebl 
12001; Volonteri & Rees 2006). Such efficient accretion would 
mean that "miniquasars" may have been abundant during 
reionization. 

Feedback from black hole accretion is a very active field, 
focusing mainly on its relationship to galaxy formation (e.g.. 
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ISilk & ReesI 119981; ICiotti & Ostrik e^ 1200 Ih. A typical ap- 
proach is that described by Springel et al T (l2005h . where some 
fraction of the rest-mass energy accreted onto the black hole 
is deposited as thermal energy within a finite volume of sur- 
rounding matter as an additional term i n the energy e qua- 
tion. This technique was recently used by lLi et all (l2007h and 
iPelupessv. Pi Matteo. & Ciardil OOOTh to study the growth of 
Pop III remnant black holes. While giving insight into the 
late time growth of the seed black holes, these simulations by 
design do not address the earliest phases of accretion. 

In this Letter we focus on the evolution of remnant black 
holes immediately following the death of their progenitor Pop 
III stars. Our simulations account for the ionizing radiation 
from the Pop III star itself, which photoevaporates the gas 
within its host halo. Because the simulations are run from 
cosmological initial conditions, we can follow the fate of the 
black hole as it encounters cold, dense gas from nearby halos, 
as well as the recombining gas in the relic H II region. We 
model the ionizing radiation from accretion onto the black 
hole by full radiative transfer coupled to the hydrodynamics, 
revealing how it regulates black hole growth and alters the 
nearby environment. We present our simulation method in §2 
and our results in §3. We conclude with a discussion in §4. 
We will present the results of two otherwise identical simula- 
tions runs, one in which ionizing radiation from the accretion 
is included ("feedback"), and one in which it is not ("no feed- 
b ack"). The simulation s with no feedback were first presented 
in lAlvarez etal] (l2008h . 

2. METHOD 

We use the code Enzo^ (lO'Shea et al.1 1200 4') to simulate 
the formation of the first star within a periodic volume 250 
comoving kpc on a side. The dark matter particle mass 
within the central ^ 90 comoving kpc, where refinement is 
allowed to proceed adaptively, is 3OM0. A fully-formed 100 

^ http://lca.ucsd.edu/projects/enzo 
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Fig. 1. — Density-squared weighted rojections of temperature (top two rows) and liydrogen number density (bottom two rows) of cubic regions 87.5 (large 
panels) and 6.25 (insets) comoving kpc on a side, centered on the position of the black hole at (left - z ~ 18) 100 (center - z ~ 14) and 200 (right - z ~ 11) 
Myr after the formation of the black hole. Top rows are without feedback, while lower rows are with feedback. While the large scale structure of the volume is 
hardly affected by the feedback, the insets, on scales of approximately 300 proper pc, show a substantial heating and density reduction within the host minihalo. 



Mo-Star is assumed to form when the overdensity is 5 x 10^ 
and the molecular hydrogen mass fraction is greater than 
5 X 10"^. The star is assigned an ionizing photon luminosity 
of 1.2 X 10^" s~' and shines for its main sequence lifetime of 
2.7 Myr Adaptive ray-tracing is used for the radiative transfer 
of ionizing radiation, which is coupled to the hydrodynamics 
via the energy equation and integrated into Enzo's primor- 
dial chemistry solve r. For more detail s on the simulation set 
up and method, see iAbel et al.l (l2007b (the set up is just the 
same as there, but with a different set of random modes for 
the initial conditions). After the lifetime of the star, radiation 



is turned off and it is assumed to collapse directly to a black 
hole with the same mass, position, and velocity as the progen- 
itor. We track the position of the black hole as the relic HII 
region evolves for 210 Myr from z~ 18toz~ 11. 

We base our model for the spectrum of accreting gas on an 
approximate form of the multi-color disk -i- power-law model 
adopted by Kuhlen & Madau (2005). We take (x u for 
hv < 200 eV, cx z^"' for 200 eV < hu < 10 keV, and L^ = 
for hi'> 10 keV. This implies a mean photoelectron energy of 
460 eV, which we take as our monochramatic photon energy. 
While admittedly crude, our model is adequate in this context. 
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Fig. 2. — Density and temperature at position of black hole along its path 
(starts at upper right) with (solid) and without (dotted) feedback. The points 
along the solid line beginning at the top right of the diagram are during the 
progenitor star's lifetime, and each point thereafter is spaced approximately 
2 Myr apart. Also shown as dashed lines are loci of fixed Bondi accretion 
rate, in units of the mass accretion rate if the black hole were shining at the 
Eddington luminosity with radiative efficiency, e = 0.1. 

where we focus on the larger scale environment on 1 0- 1 000 pc 
scales. Neglecting dissociating radiation places an upper limit 
on positive feedback effects, and is therefore conservative in 
this respect. The normalization of the total luminosity is ob- 
tained by assuming that ten per cent of the rest mass energy 
of matter falling in is converted to radiation, an d the accre- 
tion rate is determined from the Bondi-Hoyle (Bo ndi & Hovld 
fT944t lBondilfT95"2h formula, Mbh = 4TrG^M^ p/(c^, + vl^y\ 
where Cs is the sound speed of the gas and Vrei is the velocity 
of the black hole relative to nearby gas. The Bondi radius, 
tb ~ 10'^ cm (M./1OOM0)(1OOO K/T), is close to our reso- 
lution limit, indicating the local conditions at the black hole 
are appropriate for determining the accretion rate. Snapshots 
of the simulation are shown in Fig. [1] 

3. RESULTS 

The black hole is tightly bound to the center of the dark 
matter halo in which its progenitor formed for the duration 
of the simulation. Its path through density-temperature space 
is shown in Fig. |2l in cases with and without feedback from 
accretion included. Also plotted are lines of constant Bondi 
accretion rate in units of the Eddington accretion rate, Mndd = 
AirGMnipil - e)/(eaTc), where the radiative efficiency is set 
to e = O.L For M = 100 Mq, the Eddington accretion rate is 
~2x lO-^Meyr"'. 

Initially, the black hole is bathed in hot, ionized gas, with a 
temperature ^ lO'* K and hydrogen number density ^ 1 cm"^. 
The photoevaporative wind from the first star continues out- 
ward under its own inertia with a velocity ^ 10 km s"'. The 
black hole eventually encounters a nearby halo with which its 
own host halo is merging, after which the evolution depends 
on whether or not radiative feedback due to accretion is in- 
cluded. Without feedback included, the path through density- 
temperature space oscillates between cold dense gas close to 
the neighboring halo's center, to the warmer more rarefied gas 
in the ambient relic H II region, due to the orbital motion 
of the merging halos. With feedback included, however, the 
black hole stays in relatively low-density, high-temperature 
gas for the entire time, partially ionizing and heating gas 
within hundreds of pc (Figs. [T]and|3]l. 
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Fig. 3. — Profiles of hydrogen number density (top), ionized fraction (mid- 
dle), and temperature (bottom), along line connecting the black hole (at r = 0) 
to the highest density point nearby, with (right) and without (left) feedback. 
The linestyle distinguish times after the black hole is formed. Their times are 
given in the lower right panel. 

Shown in Fig. |4] are distance to the nearest high density 
peak ("clump"), the heating time at the clump, and accretion 
rate (in units of MEdd ^ 2 x 10"^ Moyr"') . This figure illlus- 
trates the in-spiralling motion of the clump and black hole as 
their host halos merge. When the clump approaches, the ac- 
cretion rate rises, accompanied by a shorter heating timescale. 
The dotted line in the top panel is the sound crossing time of 
a 1 pc clump at ~30,000 K. Before the clump approaches 
the black hole for even the first time, the heating timescale 
becomes shorter than its sound crossing time, and photoe- 
vaporation ensues. For a radiative efficiency of ten percent, 
self-regultation of accretion onto black holes within miniha- 
los occurs at accretion rates some five orders of magnitude 
below the Eddington rate, at about lO~"M0yr"'. 

4. DISCUSSION 

We have presented results from a simulation of radiative 
feedback from accretion onto a remnant black hole of an early 
Pop III star. This feedback makes it even more difficult for 
the black hole to encounter cold, dense gas, further depleting 
the reservoir of gas available for accretion. Including radia- 
tive feedback lowers peak accretion rates from ~ lO~^M0yr"' 
to ^ 10~"MQyr"', indicating that radiative feedback is very 
efficient at halting accretion. Furthermore, the feedback is 
sufficiently strong to prevent further cooling and star forma- 
tion within its host minihalo, inhibiting star formation in the 
halo over the 200 million years we have evolved it. Further 
star formation can only occur in larger systems, for which the 
dark matter potential well is deep enough or within which star 
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Fig. 4. — Evolution of heating time (top), distance from black hole to 
clump (middle), and accretion rate (bottom). Solid lines are with feedback, 
dotted lines are without feedback. The horizontal Une in the top panel is the 
sound crossing time across a 1 pc photoionized clump, while the one in the 
lower panel is the Eddington accretion rate for a 100 Mq black hole. 

formation can occur outside of the center of the halo. 

These resuhs are in agreeement with previous work, which 
indicated that massive Pop III stars forming in miniha- 
los with mass ^ IQ^Mq photoevaporate the halo gas from 
the inside out, leaving any remnant bl ack hole starved, at 
least momentar i ly, of a ccreting material jWhalen et al]|2004) : 
Kitava ma etaP l2004, i Q'Sheaet all l2005h . In particular, 
Johnson & Bromm ( 2007) showed that this photevaporation 
was likely to delay significant accretion for on the order 
of 100 million years, creating an early "bottleneck" in the 
black hole's growth. Our results imply that the situation 
is likely to be even more drastic, because radiative feed- 
back is so effective at preventing gas from reaching the cen- 
ter of the halo, even without taking into accou nt the effect 
of radiation pressure (e.g ., Haehnelt 1995; Diiks tra & Loebl 
^008; Milosavlj evic et al. 2008). An appealing alternative 
scenario for seed bl ack hole formation an d accretion was re- 
cently discussed by iDiikstra et al.l (l2008l) . in which interme- 
diate mass black holes can form by direct collapse in halos 



(e.g..l Bromm & Loebll2003l) . More detailed numerical simu- 
lations of the kind presented here will be required to test this 
hypothesis. 

The possibility that accretion onto these black holes could 
power miniquasars has received a great deal of attention 
lately, since their radiation could have significantly increased 
the ionizing photon b udget of the universe and extended the 
reionization process jRicotti & Ostrikej 120041 : iMadau et alj 
l2004l) . Because their radiation would have a long mean 
free path, they would also change the topology of reion- 
ization, blurring the distinction between ionized and neutral 
gas. They could also change the pattern of fluctuations in 
the 21 -cm background befo re reioinization, as pointe d out by 
IChuzhov et all |2006) and IChen & Miralda-Escudd (|2006|) . 
Black hole remnants of Pop III stars that formed in miniha- 
los are not likely to be the origin of the miniquasars. As we 
have shown, significant accretion onto these black holes will 
be delayed until they fall into more massive dwarf galaxies, 
which are likely to be the sites of significant star formation, 
perhaps of Pop II stars. Black hole remnants from these ob- 
jects are likely to be both more numerous, and in locations of 
higher gas density than the black holes that formed in mini- 
halos, thus replacing the Pop III seeds in models like that of 
iTanaka & HaimanI ( l2008h . 

In future work it will be necessary to determine at which 
values of the radiative efficiency feedback from black hole 
accretion becomes irrelevant. More importantly, we must bet- 
ter constrain the angular distribution (i.e. beaming), intensity, 
and spectral shape of the radiation produced by the accretion 
flow around massive early stellar mass black holes in the en- 
vironments we detailed here. Such modeling will require far 
more spatial and temporal dynamic range than has hitherto 
been feasible. 
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